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Abstract 

Context. There are only a few tens of open clusters for which ellipticities have been determined in the past. 
Aims. In this paper we derive observed and modelled shape parameters (apparent ellipticity and orientation of the 
ellipse) of 650 Galactic open clusters identified in the ASCC-2.5 catalogue. 

Methods. We provide the observed shape parameters of Galactic open clusters, computed with the help of a multi- 
component analysis. For the vast majority of clusters these parameters are determined for the first time. High resolution 
("star by star") N-body simulations are carried out with the specially developed ^GRAPE code providing models of 
clusters of different initial masses, Galactocentric distances and rotation velocities. 

Results. The comparison of models and observations of about 150 clusters reveals ellipticities of observed clusters which 
are too low (0.2 vs. 0.3), and offers the basis to find the main reason for this discrepancy. The models predict that after 
« 50 Myr clusters reach an oblate shape with an axes ratio of 1.65 : 1.35 : 1, and with the major axis tilted by an angle 
of qxY ~ 30° with respect to the Galactocentric radius due to differential rotation of the Galaxy. 

Conclusions. Unbiased estimates of cluster shape parameters requires reliable membership determination in large cluster 
areas up to 2-3 tidal radii where the density of cluster stars is considerably lower than the background. Although 
dynamically bound stars outside the tidal radius contribute insignificantly to the cluster mass, their distribution is 
essential for a correct determination of cluster shape parameters. In contrast, a restricted mass range of cluster stars 
does not play such a dramatic role, though deep surveys allow to identify more cluster members and, therefore, to 
increase the accuracy of the observed shape parameters. 
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1. Introduction lation to cluster ellipticity was briefly discussed. In the 

present study we discuss in detail all the issues related to 
Our current project aims at studying the properties of ^j^g ghape of local clusters 
the local population of Galactic open clusters. The sam- 
ple contains 650 open clusters and cluster-like associa- The shape of a star cluster, as well as its size and mass 
tions identified in the all-sky compiled c atalogue of 2.5 are the most important dynamical parameters. They are 
million stars ASCC-2.5 (|Kharchenkoll200ll ). For each clus- predesignated already at the stage of the cluster forma- 
ter, a combined spatio-ki nematic-photometric mem bership tion, when the cluster keeps the memory of the size and 
analysis was performed (jKharchenko et all l2004l ) and a shape of the parent cloud. In addition to internal pro- 
homog eneous set of diffe rent cluster parameters was de- cesses (self-gravitation and rotation) a considerable role 
rived pCharchenko et al.l 2005a. b) . In two re cent papers in cluster shaping is played by external forces, e.g. by the 



(jPiskunov et all l2007t iPiskunov et all [20083) tidal radii Gak c tic tid al field, which is acting in two ways (cf. IWielenI 



and masses of open cluster were determined and their re- 11974 Il985( l: a) stretching the cluster into an ellipsoid di- 
rected towards the Galactic centre, and b) producing clus- 
Send offprint requests to: R.-D. Scholz ter tails outpouring from the ellipsoid endpoints (cluster 

* The determined shape parameters for 650 clusters are listed Lagrangian points). The tails are composed of stars lost by 
in a table that is available in electronic form at the CDS the cluster, which lead and/or trail the cluster along its 
via a nonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via orbi t due to differential roi^ation of the Galactic disk (see 



http://cdsweb.u-strasbg.fr/ cgi-bin/qcat?J/A-|-A/| e.g. IChumak fc Rastorguevl l2006f ) . The encounters of star 



clusters with giant molecu lar clouds random ize the regular 
effect of the Galactic field (|Gieles et al.ll2006t ). Additionally, 
molecular clouds produce a screening effect, when they 
partly overlap the clusters. 

Most of the above reasons lead to the violation of clus- 
ter's sphericity with different effects on the cluster core and 
corona. The analysis of such violations and their compar- 
ison with predictions of N-body models can shed light on 
the dynamical history of open clusters. We assume in the 
following that open clusters can be represented by triaxial 
ellipsoids. 

The violations of the sphericity of globular clusters are 
evident and their ellipticities were determined for the first 
time al ready at the begin n ing of the last century. For ex- 
ample, IShaplev fc Sawverl (jl927D published the estimates 
of ellipticity of 75 globular clusters. Currently, shape pa- 
rameters ar^^etemimfidfo^ hundred Galactic globular 
clusters (jWhite fc Shawll[l987f l. i.e. according to Harris' on- 
line catalogu43 for about two thirds of the known objects 
in the Galaxy. 

The current status of the shape measurements for open 
clusters is much poorer. Until now, indications of a fiat- 
tening were obtained for a few tens o f open cluster s only . 
Raboud fc MermilliodI (|l998allbl ) and lAdams et all (|200ll 



20021) found evidence of a flattening in the Pleiades and in 



Praesepe. Their results support the predictions of Wielen's 
model that explains the ellipsoidal form of clusters by tidal 
coupling with t he Galaxy. In contrary, the flattening found 
in the Hyades (IOortlll979l) de viates from theoretical expec- 
tations. Recently, Chen et al] (2004) published morphology 
parameters, including cluster ellipticities, for 31 Galactic 
open clusters, located preferentially in the anticentre direc- 
tion of the Galaxy. 

With this paper we fill this gap in the parameter list 
of open clusters and determine the sha pe parameters for 
all 6 50 open clusters in our sample (jKharchenko et ahl 
I2005ajbl . For a comparison with these observed shape pa- 
rameters we carried out a set of high resolution ("star 
by star" ) N-body dynamical calculations of cluster models 
with different initial angular momenta located at different 
Galactocentric distances in the Milky Way. 

The paper has the following structure. In Sec. [2] we 
provide the basic equations used in the present work; in 
Sec. [3] we describe the N-body models; in Sec.|4]we compare 
the observed shape parameters with model calculations and 
with literature data; in Sec. [5] we summarize the results. 



Currently, only the two nearest clusters (Ursa Majoris 
and the Hyades) provide sufficiently accurate individual 
distances of their stars allowing for a direct determination 
of the 3D structure parameters of these clusters. For the 
other clusters we can only consider the stellar distributions 
in projection onto the celestial sphere. 

To derive the shape parameters of clusters: ellipticity, 
lengths and directions of the ellipse axes, we applied a mul- 
ticomponent analysis of the positions of cluster members in 
the standard (or tangential) coordinate system x, y in each 
sky area containing a cluster. For each cluster member, x, y 
are computed from its Galactic coordinates l,b as: 

sin(^ — Ic) 



y 



tan b ■ sin be + cos be ■ cos{l — Ic) ' 

cos be ■ tan b — sin be ■ cos{l ~ Ic) 
sin be ■ tan b + cos be ■ cos{l — Ic) ' 



where Ic, be are the Galactic coordinates of the cluster cen- 
tre. The axis x is parallel to the Galactic plane, and the axis 
y is pointing to the North Galactic Pole, positive directions 
of X, y coincide with positive directions of ^, b, and for the 
cluster centre = 0, j/c = 0. The 2nd order momenta of 
standard coordinates 



N 



N 



N 



where i — 1,2, ...N, and N is number of cluster members, 
are the basis for the characteristic equation: 



M. 



xy 



Myy-A 



0. 



The roots of this equation: A = ^/Ai, B = VA2 scale the 
principal semi-major and semi-minor axes of an apparent 
ellipse given by the distribution of cluster members over the 
sky area. An apparent ellipticity e (ellipticity hereafter) is 
then computed as 

The orientation of the ellipse is defined by the parameter 
q, which is the angle between the Galactic plane and the 
largest principal axis of the ellipse: 



tan(2(7) 



2- 



(2) 



2. Calculation of observed parameters of cluster 
shape 



and varies between 0° and 90°. 

The rms errors of these values are computed from: 



In this study we use the following coordinate sys- 
tems: spherical Galactic coordinates {l,b); the rectangu- 
lar Galactocentric system {X,Y,Z), with origin in the 
Galactic Center, and axes directed to the Sun - X, along 
the Galactic rotation at the Sun's location - Y, and to the 
North Galactic Pole - Z; and a proper rectangular coor- 
dinate system {X', ¥', Z') with origin in the centre of a 
cluster under consideration. The AT'-axis is directed along 
the projection of the cluster Galactocentric vector onto the 
Galactic plane, Y' points to the direction of Galactic rota- 
tion and Z' is directed to the North Galactic Pole. 



^ http://www.physics.mcmaster.ca/resources/globular.html 
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Figure 1. Aspect angles A and /3 versus Galactic coordi- 
nates I and b, respectively, for all 650 clusters. 



and ea, £b (as well as Eg and Sq) are calculated applying 
the law of error propagation. 

If the orientations of cluster ellipses are not random, 
one expects observing a longitudinal dependence of e. Since 
clusters reside at different heliocentric distances d, it would 
be more convenient to consider a related angle A instead 
of I, where A is the angle between the projections of the 
heliocentric and Galactocentric vectors of a cluster onto 
the Galactic plane, and can be determined from: 



cos(180° - A) = 



d — Rp^ ■ cos L 



— 2- Rp; ■ d ■ cos L 



The angle A varies between 0° and 180°. Similarly, (3 is the 
angle between the projections of the two vectors onto the 
meridional plane of a cluster: 



cos(180° - 13) 



d — Rp^ ■ cos br 



R% 



— 2- Rq ■ d ■ cos b. 



the angle f3 varies between 0° and 90°. Here d is the distance 
of the cluster from the Sun, Rq = 8.5 kpc is the distance 
of the Sun from the Galactic centre. We define A and /3 to 
be A = 180° at / = 180° and ^ = 0° at 6 = 0°. In Fig. □ we 
show the values of A and f3 for all 650 clusters. 

Using the equations above, one can compute shape pa- 
rameters (principal semi-axes A and B, ellipticity e and 
angles q, (3), both for real and modelled clusters. 



3. Numerical modeling 

A non-sphericity of star clus ters is p r edicted both by theory 
and numerical simulations. IWielenI (|l974l . Il985l ) predicted 
that the ratios of the three orthogonal axes of the cluster 
ellipsoid should show the ratio {a : b : c) = (2.0 : 1.4 : 1.0). 
The largest axis is pointing to the Galactic centre and the 
smallest one is directed to the North Pole. These results 
were obtained with N-body calculations of 500 particles, 
and were confirmed later with some what more populated 
models including up to 1000 particles (|Terle vich 1987), and 
up to 2500 particles (jChumak fc Rastorguev. ,2006 ). Such 
a number of particles corresponds to a cluster initial mass 
below IO^Mq. Presently there are evidences, however, that 
averag e masses of form i ng star clusters are larger. For ex- 
ample, iPiskunov et al.l (|2008bf ) have found, that the aver- 
age mass of the Galactic open clusters at birth is equal 
to 4.5 • 10^ m0. Therefore, it is necessary to consider more 
populated models than it has been done before. 



Table 1. The parameters of the Galactic potential compo- 
nents. 



Mass component 


M/m© 


a, kpc 


b, kpc 


Bulge 


1.4 


10^° 


0.0 


0.3 


Disk 


9.0 


10^° 


3.3 


0.3 


Halo 


7.0 


10" 


0.0 


25.0 



3.1. <j)GRAPE N-body code 

For our high resolution N-body simulations we use the spe- 
cially developed i^GRAPE code. The code itself and also 
the G RAPE hardware we used are described in the pa- 
per by iHarfst et al.l (|2007t ) in more detail. Here we just 
briefly mention a few of the special features of our code. 
The program was already thoroughly tested with different 
N-body applications, including the high resolution, direct 
study of the dynamical evolution of the Galactic centre with 
Binary (or Single) Black Holes (jBerczik et al.ll2005l [20061 : 
iMerritt et aT]l2007D Fl. 

The program acronym (/)GRAPE means: Parallel 
Hermite Integration with GRAPE. The serial and parallel 
versions of the program are written from scratch in ANSI-C 
and use the standard MPI library for communication. For 
the integration of the star cluster's dynamical evolution in- 
side the Galactic potential we use the parallel GRAPE sys- 
tems developed at ARI Heidelberg (GRACE - year 2005) 
and at MAO Kiev (GRAPE/GRID - year 2007). 

The code uses the 4-th order Hermite integration 
scheme for the particles with the hierarchical individ- 
ual block timesteps, together with the parallel usage of 
GRAPE6a cards for the hardware calculation of the accel- 
eration a and the first time derivative of the acceleration 
a (this term is usually called "jerk" in the N-body commu- 
nity). 

We specially check the effect of N-body softening, which 
suppresses the hard binary formation in our code, on the 
evolution of the cluster model. We vary the softening pa- 
rameter from a few hundreds down to a few astronomical 
units and do not find a significant difference in the evolution 
of the cluster mass and shape. We also consider the role of 
a possible initial mass segregation and do not observe any 
dependence of the evolution of the semimajor axis on this 
value. 

Compared with the previous public version we add two 
major changes to the code. First of all we add the pos- 
sibility to have some external potential, acceleration and 
also "jerk" . For the external po t ential we choose the form 
proposed bv lMiyamoto fc Nagail (|l975f ). Using such a multi- 
mass component potential we can easily approximate the 
Galaxy's external force, acting on our star cluster in the 
Galactic disk at different Galactocentric distances. The sec- 
ond change includes the possibility to turn on the mass loss 
due to the stellar evolution for every modelled star parti- 
cle. For mctallicity-dependent stellar lifetimes we use the 



^ The present version of the code and the full snapshot 
datasets analyzed in the paper are publicly available from: 
ftp : //ftp . ari .uni-heidelberg.de/pub/ staff/ 
berczik/phi-GRAPE-cluster/. 



Table 2. The list of model parameters: initial mass Mc(0), 
number of particles N, cluster radius Rc, distance from 
the Galactic center Rq, circular velocity Vb, time step St, 
dimensionless angular velocity ujq; concentration parameter 



Wo = 6.0. 


Mc(0), 


iV 




-Ro, 


Vb, 


5t, 




rriQ 




pc 


kpc 


km/s 


Myr 




10^ 


4040 


3 


7.0 


236 


2.45 


0.0, 0.3, 0.6 


10^ 


4040 


3 


8.5 


233 


2.45 


0.0, 0.3, 0.6 


10^ 


4040 


3 


10.0 


231 


2.45 


0.0, 0.3, 0.6 


5- 10^ 


20202 


7 


7.0 


236 


3.91 


0.0, 0.3, 0.6 


5- 10^ 


20202 


7 


8.5 


233 


3.91 


0.0, 0.3, 0.6 


5- 10^ 


20202 


7 


10.0 


231 


3.91 


0.0, 0.3, 0.6 


10^ 


40404 


10 


7.0 


236 


4.71 


0.0, 0.3, 0.6 


10'' 


40404 


10 


8.5 


233 


4.71 


0.0, 0.3, 0.6 


lo-* 


40404 


10 


10.0 


231 


4.71 


0.0, 0.3, 0.6 



approximation formula proposed bv lRaiteri et al.l (|l996l ): 

logT = ao{Z) - ai{Z) • logm + a2{Z) ■ log^ m, 

where r is expressed in years, the stellar mass m in solar 
masses, and where Z is the abundance of heavy elements. 
The coefficients are defined as: 

log' Z, 



aoiZ) = 10.130 + 0.0755 -logZ- 0.0081 
ai {Z) ^ 4.4240 + 0.7939 • logZ + 0.1187 • log^ Z. 
a-iiZ) = 1.2620 + 0.3385- log Z + 0.0542 



log' Z. 



The mass lost by stars during their evolution is cal- 
culated with the help of tables from the paper of 
Ivan den Hoek fc GroenewegenI (|1997[ ) and is approximated 
analytically via the metallicity-dependent formula: 



''eject 



-0.4205 -Z 



-0.0177 



(0.9015 + 0.6294 -Z) ■ m,. 



For simplicity, we assume in our model that stars loose their 
masses permanently with a constant rate: —nif.ject/T. 

Since stars lose the bulk of their masses at the end 
of their evolution, we test, also, an alternative mass loss 
scenario, where we reduce the stellar masses just in one 
timestep at the end of the star life. We find that the mod- 
els show almost the same evolutionary patterns for shape 
parameters and mass loss of clusters independent of the 
version of stellar mass loss used. 

3.2. Initial conditions for the star cluster 

For the generation of the initial position and velocity distri- 
butions of cluster particles we use the standard King model. 
The code for the generation of the initial rotating star clus- 
ter is described in detail bv lEinsel fc SpurzemI (|l999f l. We 
also use their notation to parametrize our rotating King 
model family with the concentration parameter Wo and the 
dimensionless angular velocity ujq. Thus each model can be 
parametrized by the combination of these two numbers. 
The larger Wo the larger is the central concentration. The 
larger luq the larger is the angular momentum of the model. 

For the present simulations we use the variation with 
three basic King models only. We vary only the rotation 
parameter ujo = (0.0, .3, 0.6). According to Table 1 in 
lEinsel fc SpurzemI (jl999t ) the ratio of the "pure" rotational 
energy to the total kinetic energy of the model for these 



three cases is Erot/Ektn = (0.0, 0.105, 0.278). The con- 
centration parameter was always set to Wo = 6.0 (which 
corresponds to the models wit h medium conc entration). 

As a next step, taking the ISalpeteil (|l955f ) initial mass 
function (IMF) /(m) — dN/dm, we create the random par- 
ticle mass distribution: 

dn(rn) — f{m) ■ dm — C ■ m^'-^^"-' • dm, 

with the lower and upper mass limits m; — 0.08 m^, 
mu — 8.0 toq and with Salpeter's slope a = 1.35. We 
use the upper mass limit of 8.0 m©, because we consider 
a "pure/classical" open cluster (older than 30. ..40 Myr), 
when all high mass stars have already finished their life as 
a SNII, and swept out all the residual gas left from the 
cluster formation process. 

The necessary number of particles for a given initial 
mass of the cluster Mc(0), can be computed from the as- 
sumed IMF as 



N = M,(0) 



1 



(m, 



With the adopted IMF parameters for the model cluster 
and initial masses Afc(O) = (10'*, 5 • 10'^, 10'^) m© we have 
the following particle numbers N = (40404, 20202, 4040). 

After our setting of initial masses, positions and veloc- 
ities of every particle in the model, we do the standard 
N-body normalization (jAarseth et al.l 0*9741 ) of the model, 
and rescale the initial cluster data to put our system into 
virial equilibrium Egra — —'^■Ekin = '^■Etot with the follow- 
ing parameters: G = Mc = 1 and Etot = —0.25. This step 
does not change any physical quantities of the modelled star 
cluster, but is just very useful for numerical reasons. For a 
King model with Wq — 6.0 such a normalization produces 
a half-mass radius of about 0.8 (in dimensionless units). 

After the construction of dimensionless parameters of 
our model we can easily extend them to a different phys- 
ical cluster mass and half-mass radius. In principle, these 
two parameters can be set independently, but in order to 
reduce the number of independent initial parameters we 
decide to use some physically motivated relation between 
initial mass and initial half-mass radius. For this purpose 
we use the extension of the well known mass vs. radius re- 
lation for molecular clouds and clumps in the Galaxy (see 
iTheis fc Hensleilll993t llnoue fc Kamavall2000( l. which can 
be scaled as: 



Rc ~ 100 ■ 



Mr 



10*3 rriQ 



pc 



For the three physical masses chosen we set the correspond- 
ing radii to Rc — (10, 7, 3) pc. 

In order to check the dependence of the computed evo- 
lution of cluster shape on the adopted mass vs. radius rela- 
tion, we run a series of models where we set the normalising 
factor to 100, 71, and 43 pc. We find that for typical model 
parameters, the evolution of the cluster shape does not de- 
pend on this factor, in practice. 

3.3. Galactic rotation curve 

For the external Galactic potential shaping a cluster we 
choose the combined "Plummer-Kuzmin disk" form (see 



1000 me 



o 



5000 mp 







10000 mz 



































































i 




















































































































































































- 






"0= 


0.0 














0.0 














0.0 




- 


































































































V 


>^ 












































- ^ 










v\ 

\x 


















\v 
















































































R 


3=8 


5 1< 




\ 








3=8 


5 k 


pc 








R 


3=8 


5 1< 


pc 














1 ^ 































0.0 0.2 0.4 0.6 O.B 1.0 0.0 0.2 0.4 0.6 0.6 1.0 0.0 0.2 0.4 0.6 0.8 1.0 

t, Gyr t, Gyr t, Gyr 

Figure 2. The evolution of the mass of dynamically bound 
particles for three cases of the initial mass of the cluster 
models (three columns with masses shown at the top of each 
column) . The upper panels are constructed for non-rotating 
models {uq = 0.0) and show the effect of the Galactocentric 
distance. The dotted (red), solid (black) and dashed-dotted 
(blue) curves are models at Rq ~ 7, 8.5 and 10 kpc. The 
bottom panels illustrate the effect of cluster rotation at the 
Solar Galactocentric distance. The solid (black), dashed 
(green) and long dashed (magenta) curves show models 
with Wo = 0.0, 0.3 and 0.6. 



iMivamoto fc NagailfTOTSl l: 
<i>(r, z) 



G-M 



^Jr^ + {a + ^J\P- + 

Coupling such a potential with a three-component Galactic 
mass distribution model comprised of "Bu lge", "Disk", 
"Halo" components llDoupho le fc Colinlll995 ^ one can eas- 
ily reproduce the observed rotation curve of the Galaxy 



1/2 (r) 



-V$B«i3e(r, Z) - \/'!>Dtskir, Z) - V^Halo{r, z) . 



In order to bring this into agreement with the kinemat- 
ics of the open cluster subsyste m in the Solar neighbor- 
hood we slightly modify the Doup hole fc ColinI (jl995f ) pa- 
rameters. The parameter values consistent with the ob- 
served Oort constant s {A,B) = (14. 5 ± 0.8 , -13.0 ± 1.1) 
km/s/kpc derived by iPiskunov et aT] (|2006f ) are shown in 
Table. [TJ The adopted values of Oort constants correspond 
to V/Rq ^ 27.5 ± 1.3 km/s/kpc, which at Rq = 8.5 kpc 
gives V w 233 km/s. Using these data we derive circular 
velocities at distances of 7.5 kpc and 10.0 kpc from the 
Galactic centre, too. 



3.4. Model list and model cluster memberships 

As a starting point for our model cluster, we select the 
position inside the Galactic disk {Rq, 0.0, 0.0) with the 
corresponding circular velocity (0.0, —Vq, 0.0). For the set 
of our runs we use three values for the positions inside the 
Galactic disk: Rq = (7.0, 8.5, 10.0) kpc with correspond- 
ing circular velocities Vq = (236, 233, 231) kni/s. Together 
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Figures. The particle distribution in the X'Y'— and 
Z'Y'— planes of a model with initial mass Mc{0) — lO'^m© 
and distance from the Galactic centre Rq = 8.5 kpc at dif- 
ferent moments of time. The cluster age is indicated in the 
right corner of each panel. Cyan dots - all particles, red 
dots and red-color labels - dynamically bound particles, 
black dots and labels - massive particles (m > Im©), black 
labels in brackets - dynamically bound massive particles. 
Arrows indicate the limits of cluster tails when they do not 
fit to the frame. 



with the three initial masses Mc(0) of clusters and with the 
three rotational parameters wq these three positions create 
the set of 3 X 3 X 3, i.e. in total 27 models (see the full list of 
model parameters in Table [2|). Each model comprises a few 
hundreds of "snapshots" made at different moments of time 
equally spaced from t — Otot > 1 Gyr. For each model par- 
ticle, every snapshot displays its initial and current mass, 
coordinates and velocities in the rectangular Galactocentric 
coordinate system. 

To provide the detailed analysis of star cluster shapes 
we need to find the cluster centre position, at first. This 
issue requires special attention, since particles having left 
the cluster during the evolution should not be taken into 
account for this purpose. For this we use our own iterative 
routine considering the distributions of positions and ve- 
locities of particles and apply it to every snapshot. We find 
that we are able to determine the centre even for models 
with highest mass loss, which lose more than half of their 
initial content during the evolution. 

As a next step in our analysis, we convert the positions 
and velocities of the particles into the proper coordinate 
system of the cluster related to the local density centre 
and define the list of particles which are still bound to the 
cluster. We use the particle kinetic and gravitational energy 
as criteria of belonging to the cluster. The kinetic energy is 
computed from the velocity of every particle in the proper 
system of the cluster. For the gravitational energy we use 
the value of the cluster's self-gravity potential, determined 
as the sum of all interactions between a selected particle and 
the rest of the particles. We assume, that only the particles 
which have a negative relative energy are still bound to the 
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Figure 4. The evolution of ratios of the model ellipsoidal 
axes a : c and 6 : c for bound particles for models of different 
initial mass (three columns with masses shown at the panel 
top). In every panel the upper triplet of the curves indicates 
the a : c ratio and the lower triplet corresponds to & : c. The 
designations are the same as in Fig. [2l 

cluster: 

|£;f™| > sf". 

In other words, for bound particles we can write: 

> (3) 

Using such a condition for the definition of cluster mem- 
bership, we exclude particles which have relatively large 
velocities compared to the cluster centre and construct a 
list of particles which we select to be "dynamical" cluster 
members. 

We define the current mass of the cluster Mdt) as the 
sum of masses of all bound particles. In Fig. [2] we show 
the evolution of Mc{t) for a set of selected models. The 
evolution shows a simi lar behaviour as in the models of 
lErnst et al.l (|2007[ ) and lKim et all (|2008D . Fig. [H indicates, 
that a star cluster looses from 35% to 70% of its initial mass 
during the first 1 Gyr and the mass loss is decreasing with 
increasing distance of clusters from the Galactic centre and 
with increasing cluster rotation velocity. 

We double check our mass loss sequences by comparing 
with the results which we obtain with other fully indepen- 
dent and also publicly available N-body codes (N-body4 
(with GRAPE6) & N-body6-h-h). All the models show the 
very same evolution patterns of mass loss. 

4. Shape parameters of observed and modelled 
clusters 

4.1. Shape parameters of the cluster models 

Let us first consider how the cluster model, initialized as 
a spheroid, changes its shape and orientation with time. 
The evolution of a selected model in the planes X'Y' and 
Y'Z' is shown in Fig. [31 With time the cluster elongates 
along the line of the Galactocentric radius, and begins to 
"leak" losing low-mass particles from the limits correspond- 
ing to the Lagrangian points, farthest and nearest to the 
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Figures. The evolution of the tilt angle qxY computed 
for dynamically bound particles for three cases of the ini- 
tial mass of the cluster models (three columns with masses 
shown in the top of the panels) . The designations are the 
same as in Fig. [5J 



Galactic centre. Because of the differential rotation the par- 
ticles overtake the cluster at the end nearest to the Galactic 
centre and lag it at the farthest end. The tails contain stars, 
which are not bound to the cluster anymore. The dynam- 
ical members form an ellipsoid projected onto the X'Y' 
and Z'Y' planes as ellipses. In the X'Y' plane the ellipse is 
tilted with respect to the X'-axis at an angle qxY, whereas 
in the Z'Y'- plane it follows the F'-axis. 

Such details in the distribution of N-body particles arise 
due to interaction with the Galactic tidal field. They are 
found in all N-body simulations of open clu sters (see, e.g. 
Fig.7 of lTerlevichlll987t IChumak fc Rastorguev.20061 . 

In Fig. [3] one can also clearly see the so-called "tidal tail 
clumps" (star density enhancements) at about 150 pc from 
the cluster centre. Such clumps were men t ioned probably 
for the first time in lCapuzzo Dolcetta et aLl(l2005l) and wer e 
discussed in more detail recentlv bv lKiipper et al.l (|2008h . 
We observe such clumps in all our models (see the model 
videos at the FTP link shown in the foo tnote on page[^. 
They can be explained by a simple theory (|Just et al. 200a). 

In Fig. m and Fig. \E\ we show the evolution of the axis 
ratios a : c and b : c and of the tilt angle qxY for selected 
models. Already after a few tens of Myrs the initially spheri- 
cal cluster transforms into an ellipsoid tilted with respect to 
the Galactocentric radius, with the semi-major axis about 
twice as large as the minor one, and a tilt between 30° 
and 40°. The shape parameters are almost independent of 
the Galactocentric distance of the model clusters, but the 
ellipsoid becomes fiatter with increasing cluster rotation. 

N-body models offer good possibilities for studying var- 
ious biases which can occur in cluster ellipticities derived 
from observations. In Fig. [S] we show the dependence of the 
ellipticity of a model cluster on sampling constraints due to 
the mass and the distances from the cluster centre of the 
particles considered. In this way we mimic the usual selec- 
tion effects which impact the observations of real clusters. 
Cutting the model particles at a certain mass limit simu- 
lates a possible bias due to a limited deepness of a survey. 
The restriction of the models by distances of particles from 
the cluster centre simulates the bias arising due to underes- 
timation of cluster size. Hereafter, we call the latter effect a 
"restricted area bias" . We consider a cluster model, where 
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Figure 6. Dependence of the shape of a model cluster on 
the sampling constraints. Bottom: cluster ellipticity e and 
semiaxes A and B as functions of lower mass limit (left) 
and radius of the considered area (right). The circles show 
the major semiaxes, the triangles are for the minor one, and 
short (red) horizontal bars mark the ellipticity. Vertical bars 
indicate the statistical uncertainties. The upper row shows 
the numbers of particles participating in the construction 
of the corresponding values. 

both major and minor semi-axes of the ellipsoid are seen 
under a right angle, i.e. the projected ellipticity is at its 
maximum. The model cluster is 150 Myr old, its current 
mass is 4.2 • 10'^ wq, and, as a consequence, its tidal radius 
is rt = 22.5. The number of dynamically bound particles is 
Nc = 19371. 

The tidal radius of the model cluster was computed from 
the current total ma ss of membe r particles with the well 
known relation from iKingI (|1962| ). This places the model 
tidal radii into the same system as the observed tidal radii 
of open clusters. Note that for about 1/3 of our clusters, 
we determined tidal radii by a fit of t he cluster radial den - 
sity distribution to a King profile (jPiskunov et all [20071 ). 
and for the remaining cluste rs we used a proper cal ibration 
to estimate their tidal radii (Pisk unov et al.ll2008 a^. Then, 
assuming that clusters fill up their potential well, the clus - 
ter "tidal" masses were computed from the iKinel (|1962[ ) 
formula for all 650 clusters. 

From Fig. E] we conclude that excluding low-mass parti- 
cles from the consideration has a much smaller impact onto 
the resulting ellipticity, than excluding the external parts 
of the model cluster from consideration. Indeed, a removing 
of low-mass particles causes the axes A and B to slightly 
decrease. Since they change, however, in coordination to 
each other, the ellipticity is about constant within the sta- 
tistical uncertainties. In contrast, with including more and 
more distant particles in the ellipticity calculation, the el- 
lipticity increases steadily. At a distance of two tidal radii, 
the corresponding ellipticity is close to the limiting value 
which is achieved if all bound particles within about three 
tidal radii are taken into account. 

This behaviour can be easily understood if one consid- 
ers radial dependences of the axes A and B. The minor axis 
B reaches its maximum value already within rt and does 
not increase anymore, whereas the major axis A continues 
to rise beyond the tidal radius. One should note that the 
cluster flatness is produced by a relatively small number of 
cluster particles: only 2% of bound particles occupy a zone 
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Figure 7. Dependence of the apparent ellipticity e on the 
probability threshold for cluster membership for the open 
cluster Stock 2. The bars correspond to mean errors. Large 
open circles indicate the apparent ellipticities based on 
the sample of stars with membership probability larger 
than l%(the left circle), larger than 14%(the middle cir- 
cle), larger than 61%(the right circle). 

beyond 1 rt, and only 0.1% of them can be found at r > 2 rt. 
In contrast, a few tens of the most massive particles, which 
are not constrained spatially, produce about the same ellip- 
ticity, as the rest of low mass members. A consequence of 
a smaller sample is, however, a considerably higher statis- 
tical uncertainty in the resulting ellipticity. A comparison 
with older models shows that the above conclusions are still 
valid for cluster ages up to one Gyr. 

In summary we conclude that the determination of clus- 
ter ellipticity is not influenced strongly by the survey's 
deepness, but rather by the sizes of the surveyed areas 
around the clusters. Restricting to areas smaller than the 
corresponding tidal radius one takes a risk to introduce a 
"restricted area bias" in the resulting cluster ellipticities, 
hence making clusters more circular than they are in real- 
ity. 

4.2. Determination of the uniform shape parameters of the 
observed and modelled clusters 

In order to compare the observations with models, one 
needs parameters adequately describing both entities. From 
observations, we compute apparent ellipticities e and ori- 
entation parameters q for each cluster of our sample using 
the approach presented i n Sec. [H and adopting the appar- 
ent cluster radius r2 from lKharchenko et al.l (j2005ai rb[). The 
observed shape parameters e and q for each of the 650 open 
clusters are listed in a table that is available in electronic 
form only. It can be retrieved from the CDS. To compute 
the parameters, we consider only the most probable mem- 
bers, with kin ematic and photometric probabilities higher 
than 61% (see iKharchenko et al.|[200^ . This decision was 
simply guided by the fact that a sample was stronger con- 
taminated if we included stars with lower membership prob- 
ability. Since field stars tend to a random distribution, the 
contamination leads to a bias: a fictitious decrease of the 
observed ellipticity. 

For illustration, we carried out the ellipticity calcula- 
tion considering cluster stars of different membership prob- 
ability. For our basic sample of 152 clusters (see Sec. 14.31 
for definition), we obtain an average ellipticity of e = 0.18 
when we consider stars with a membership probability bet- 
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Figures. Distributions of ellipticities e (left panels) and 
orientation angles q (right panels) for 152 clusters. The pan- 
els are marked with the identifiers of the set. 



ter than 61%, e = 0.14 with membership probability better 
than 14%, e — 0.12 with membership probability better 
than 1%, and e = 0.09 if we consider all stars projected on 
the cluster area. A probability threshold of 61% is a com- 
promise between a possible bias due to backgroud contam- 
ination and the number of stars included in the ellipticity 
calculation. 

On the other hand, one must make certain that stars 
with membership probability larger than 61% describe ad- 
equately the cluster properties. In Fig. [7] we show the varia- 
tions of the apparent ellipticity with the probability thresh- 
old in the case of the cluster Stock 2, for example. This 
cluster is a well populated cluster, so we get satisfactory 
statistics even for the smallest star samples containing the 
most probable cluster members. The contamination effect 
by field stars is strong if the probability threshold P is less 
than 30%, but it becomes almost insignificant (although 
still systematic) at 40% ^ P < 60%. If we consider stars 
with membership probabilities larger than 60%, the vari- 
ations in the resulting apparent ellipticity have a random 
character. This indicates that - although the sample of stars 
with membership probability over the 61% threshold may 
contain a number of field stars - they do not introduce a 
bias in the apparent ellipticity determination. A further in- 
crease of the probability threshold does not improve the 
results significantly, but the mean error of the apparent el- 
lipticities becomes larger due to poor statistics. Therefore, 
we consider the choice of a 61% threshold to be the opti- 
mum. 

From the models, we obtain the axes a, 6, c and the tilt 
angle qxY describing the shape and orientation of a cluster 
in the proper coordinate system {X' ,Y' , Z'). To compare 
them with the observed parameters e and q, we must "view" 
the models under the same conditions which are valid for 
the observed clusters. First, one has to remember that the 



limiting magnitude of the ASCC-2.5 is about V — 12, and 
even in the nearest clusters we can, therefore, observe only 
stars more massive than 0.7 Mq, i.e. one order of magnitude 
higher than the lower mass limit of model particles. Further, 
the resulting shape parameters are strongly dependent on 
the size of the area around the cluster centre included in the 
computations. Therefore, a model cluster must be scaled to 
the apparent cluster radius r2 of its observed counterpart. 

For each observed cluster, we first selected its theoreti- 
cal counterpart from Table [2] at the Galactocentric distance 
which best fits the location of the observed cluster. Then, 
for a given model we selected a snapshot with an age clos- 
est to the observed age. This snapshot was placed at the 
position of the real cluster, and the model particles were 
projected onto the face-on plane (i.e. for every particle the 
standard coordinates x,y were computed). From the com- 
plete list of model particles, we selected the actual members 
of the cluster by applying eq. ^ . We used this list to com- 
pute "true" model parameters of the cluster's shape with 
the formulae from Sec. [2] 

The models contain particles down to relatively low 
masses which can not always be detected in real observa- 
tions. In order to reproduce the observed conditions (i.e. 
to obtain "observed" model parameters) , we further dimin- 
ished the list of dynamical members selecting particles in 
the observed mass range and obeying the condition that 
the number of modelled and observed members within the 
cluster area r2 is approximately the same. 

As a result, we consider four different sets of input 
data. The first set includes real observations of clusters. 
We call this set "OC" . The following three sets are related 
to the models. The set "MC" represents the models as such, 
and reproduces "true" clusters containing all the dynami- 
cal members in the complete range of particle masses. The 
third set "MCm" includes only "massive" dynamical mem- 
bers of given clusters which are above the limiting mag- 
nitude of ASCC-2.5. Finally, the set "MCmr" completely 
reproduces the observed sample: it includes the same num- 
ber of particles selected in the same mass range Am as we 
do observe in real clusters. Similar to the observed clusters, 
the MCjnr models contain only those particles, which reside 
within the apparent cluster radius r2. Since the sets MCm 
and MCmr bridge the sets MC and OC, they are introduced 
to investigate the biases in the observed shape parameters 
which may occur due to the limitations of mass ranges and 
spatial distribution of cluster stars described in Sec. 14.11 

4.3. Comparison of shape parameters of the modelled and 
observed clusters 

For comparison with the model data we consider only clus- 
ters with more than 20 most probable members. We also 
exclude the two very populated clusters NGC 869 {h Per) 
and NGC 884 (x Per), as they are overlapping in the pro- 
jection onto the sky, and have a large number (more than 
50%) of members in common, making an accurate determi- 
nation of their individual shape parameters rather difficult. 
Hereafter, we refer to the resulting sample as the "basic 
sample" comprising 152 objects. 

According to the mass estimates of open clusters given 
in iPiskunov et al.l (|2008al ). the average tidal mass of clus- 
ters in the basic sample turns out to be about lO^m0 which 
is compar able to the average clus ter mass of 700 m© deter- 
mined in iPiskunov et al.l ( 2008b( l for the complete cluster 



Table 3. Projected semi-axes A and B in selected directions of the sky for the model of preferentially oriented cluster 
ellipsoids. 
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Figure 9. EUipticities e of 144 low- latitude open clusters 
{P < 20°) versus age t. The panels are marked with 
the identifiers of the set. The dotted vertical line marks 
t = 50 Myr. In each panel, the horizontal lines show the 
corresponding average ellipticity for clusters older than 
50 Myr. The dashed lines indicate their standard devia- 
tions. Vertical bars for a cluster at t — 2.46 Gyr show typ- 
ical rms errors of the ellipticity estimates. 



sample of 650 clusters. Therefore, we assume that the basic 
sample represents sufficiently well the typical cluster pop- 
ulation in the Solar vicinity, havi ng typical initia l masse s 
about 4.5- 10^ m© as determined in lPiskunov et al.l (|2008tJ) . 
For this reason, we select the cluster models with initial 
mass of 5 • lO'^m© as the most suitable ones for comparison 
with observations. According to Fig. HI the shape param- 
eters do not change considerably for models with larger 
initial masses, whereas the ellipticity becomes more promi- 
nent in clusters with initial masses considerably lower than 
the adopted STO'^toq. In the following we consider the non- 



rotating models of open clusters keeping in mind that with 
increasing rotation velocity clusters become also flatter. 

The distributions of ellipticities e and orientation angles 
q of the observed and modelled clusters of the basic sample 
are shown in Fig. [51 The distribution of ellipticity shows 
a clear maximum for the complete cluster models (the set 
MC) at e w 0.25 . . . 0.35, and the corresponding ellipsoids 
are elongated parallel to the Galactic plane (the orienta- 
tion angle g « 0°). For the models MCm and MCmr which 
take into account the restrictions set by observations, the 
ellipticity decreases, while the spread in q increases indi- 
cating a more random orientation of the apparent ellipses. 
The observed clusters (the set OC) show relatively small 
ellipticities with a peak between 0.075 . . .0.175. Although 
a peak at g = is still observed, the orientation angles 
q are distributed over the whole range. This is simply a 
consequence of small ellipticities since apparently circular 
clusters have no orientation angle. 

We show the dependence of ellipticities e on cluster age 
for all the four sets in Fig. El We limit our consideration to 
144 clusters with (3 < 20° in order to minimize the impact 
of the projection onto the smallest axis of the ellipsoid oc- 
curring for clusters at large galactic latitudes. According to 
the model set MC, clusters change their initially circular 
shapes into ellipsoids during the first 50 Myr and keep this 
form thereafter. This behaviour can be also concluded from 
Fig. |3l The small scatter around the average ellipticity for 
clusters older than 50 Myr is defined by their location in 
the Galactic plane relative to the Sun and to the Galactic 
centre, i.e. by the aspect angle A. When we reduce the num- 
ber of observable members by excluding model particles of 
low masses (the set MCm), the average ellipticity does not 
change, although the distribution shows a larger scatter 
which arises due to poorer statistics of the remaining dy- 
namical members. When further excluding model particles 
located outside the apparent radius assumed for the clus- 
ters, the average ellipticity becomes significantly smaller. 
The observed clusters show a similar ellipticity distribution 
as the model set MCmr • These features confirm the analysis 
carried out in Sec. 14.11 

The following statistics supports the decisive role of spa- 
tial limitation (r < r2) hiding the real ellipticity of clusters. 
The average ellipticities of 103 clusters with t > 50 Myr 
are equal to 0.328 ± 0.004, 0.322 ± 0.013, 0.184 ± 0.010, 
0.186 ±0.010 for MC, MCm, MCmr and CO, respectively 
The similar average ellipticities of the sets MCmr and OC 
indicate that a possible bias due to background contami- 
nation is rather small if the ellipticity calculation is based 
on cluster stars with membership probabilities higher than 
61%. We conclude, also, that a spatial limitation of clus- 
ter members introduces a significant bias in the determi- 
nation of ellipticities, whereas a mass limitation increases 



mainly the random uncertainties of the results. The smaller 
elHpticities of the sets MCmr and OC are results of too 
small radii r2 assumed for the clusters. Note that, on av- 
erage, the ASCC-2.5-bas ed estimates of cluster radii r2 
(jKharchenko et aTl |20053) are already larger by a factor 
of two compared to previously published values. Better es- 
timates of cluster radii are only possible if a proper separa- 
tion of cluster members from the numerous field stars can 
be achieved in outer cluster regions. This would require 
more accurate surveys of proper motions and photometric 
data all over the sky. On the other hand, surveys deeper 
than the ASCC-2.5 would increase the random accuracy of 
the ellipticity determination of open clusters. 

4.4. Distribution of the cluster ellipticity over tlie sky 

Up to now, the elliptici ty of open clusters was studied either 
for individual objects (jRaboud fc MermilliodlllQQSallbl'). or 
in se lected directions of the Galactic disk ( Chen et al.l 
l2004l) . In the present study we determined the apparent 
shape parameters for clusters observed all over the sky. 
Therefore, it seems to be useful to derive geometric rela- 
tions between the apparent ellipticity of clusters and their 
spatial location in the Galaxy. Of course, a regular depen- 
dence can only be expected if cluster ellipsoids are not ori- 
entated randomly in space but show a certain orientation 
with respect to the Galactic centre. 

Indeed, the model calculations indicate that the 
Galactic tidal field and differential rotation quickly align 
clusters along a preferential direction, and the resulting el- 
lipsoids show always the same axes ratios. Let a, 6, c be 
the semi-axes of these three-axial ellipsoids with the semi- 
major axis a tilted with respect to the Galactocentric radii 
of the clusters by an angle qxY- For simplicity, we further 
assume that the major axis is parallel to the Galactic plane. 
Hereafter, we call this approach the model of preferentially- 
oriented ellipsoids, and the ellipsoid itself the "reference" 
one. 

From observations, however, we can only see a two- 
dimensional projection of the reference ellipsoid on the sky 
that is an ellipse with the semi-axes A and B. The ratio 
B/A varies from c/a to 1, depending on the aspect an- 
gles A and /3 describing the cluster location with respect to 
the Sun and the Galactic centre. The apparent ellipticity e 
changes from to 1 — c/a, respectively. 

When a cluster is located in the Galactic plane {b = 
P — 0°), B — c is always valid, and A varies from b to a. 
The apparent ellipticity reaches the maximum e — 1 — c/a 
at aspect angles Xac = 90° + qxY (for / = 0° . . . 180°) and 
at Xac = 90° - qxY (for I = 180° . . . 360°). In contrary, the 
apparent ellipticity reaches the minimum e = 1 — c/6 at 
Abe = qxY (for / = 0° . . . 180°) and Xbc ^ 180° - qxY (for 
/ = 180° . . . 360°). For all other directions, the apparent 
ellipticities can be computed from the equations given in 
the first row of Table \3\ 

Though the majority of star clusters is located near the 
Galactic plane, there is a number of clusters at higher galac- 
tic latitudes. In a special case, when clusters are located di- 
rectly at the Galactic Poles, one observes A = a and B = b. 
For all other clusters outside the Galactic plane, the pro- 
jection of the reference ellipsoid is defined by both aspect 
angles. For example, at A = Xac the ellipticity decreases 
monotonically from 1 — c/a at /3 = 0° to 1 — &/a at /3 = 90°. 
The corresponding expression is given in the second row of 
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Figure 10. Ellipticity e versus aspect angle A for 103 clus- 
ters with (3 < 20° and t > 50 Myr. The left panel corre- 
sponds to the Galactic longitude range I = 0°...180°, the 
right panel to I = 180°. ..360°. The curves show the e(A, (3)- 
relations for /3 = 0° (upper curve) and /3 — 20° (lower 
curve) constructed for the reference ellipsoid tilted to the 
direction to the Galactic centre by angle qxY = 30°, having 
axes ratios (1.65 : 1.35 : 1). Short horizontal bars (red): set 
MC (complete models), open circles: set OC (observations). 
The vertical bars drawn for one object show the typical rms 
error of the observed ellipticity. The dashed vertical lines 
mark Xac, the dotted ones Abe- 



Table 131 At A = Xbc this relation is rather different since the 
projected ellipse changes its orientation at some angle (3o. 
Thus the ellipticity first decreases from 1 — c/b at /3 = 0° to 
at /3 = /3o, and then increases to 1 — 6/aat /3 = 90°. The 
corresponding equations are listed in the last two rows of 
Table [3l The aspect angle /3q where the projection becomes 
a circle {A = B) is determined by the semi-axes a, 6, c of 
the reference ellipsoid from the equality condition 



\/ ■ cos^ /3o + c^ • sin'^ (3q 



(4) 



For given parameters of the reference ellipsoid, the 
Galactic meridian at Xac defines a locus of maximum ellip- 
ticities e(AQc, P) over the sky, whereas the Galactic meridian 
at Xbc gives a locus of minimum ellipticities e(Ahc, P)- 

In Fig. [TO] and Fig. [11] we plot the analytical relations 
from Table [3] together with the model and observed data 
on ellipticities versus the aspect angles for "relaxed" clus- 
ters which are older than 50 Myr. Though a few observed 
ellipticities fit the predicted values, the majority of the 
observed clusters show, as expected, too low ellipticities 
due to the restricted area bias described above. Therefore, 
we cannot use the observations for deriving the param- 
eters of the reference ellipsoid. Instead, the parameters 
a, 6, c and qxY were found by fitting the relations from 
Table |3] to the data points of the complete models (MC) 
as (a : & : c) = (1.65 : 1.35 : 1) and qxY = 30°. These 
parameters are based on the non-rotating models assuming 
the same initial masses (5- IQ^Mq) but different location of 
clusters from the Galactic centre (see also Fig.|4]and Fig. [5]). 
For this reference ellipsoid, the angle /3o is determined from 
eq. gj) to be /3o = 57.6°. 
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Figure 11. EUipticity e versus aspect angle /3 for 110 clus- 
ters with t > 50 Myr. The curves show the e(A, /3)-relations 
for A = Xac (dashed curve) and A = Xbc (dotted curve) con- 
structed for the reference ellipsoid tilted to the direction of 
the Galactic centre by an angle qxY = 30°, and with axes 
ratios (1.65 : 1.35 : 1). Short horizontal bars (red): MC 
(complete models), open circles: OC (observations). The 
vertical bars drawn for one object show the typical rms 
error of the observed ellipticity. The dashed vertical line 
marks /3 = 20°. 



4.5. Comparison with the literature 

In the literature, there are o nly a few studies on the elliptic- 
ity of Galactic open cl usters. Rabo ud fc Mermilliodl (|l998a[ ) 
and lRaboud fc Mermilliod (1998b ) studied the shape of the 
Pleiades and Praesepe using the best data on astrometric 
and photometric membersh ip then avai lable for these open 
clusters. For the Pleiades. iRaboud fc M ermilliod (1998a} 
determined e — 0.17 ± 0.05. Since iRaboud fc Mermilliodl 
(|l998h ) do not provide the ellipticity for Praesepe, we com- 
puted this value (e = 0.05 ± 0.07) based on the data on the 
cluster membership from their Table 2. IChen et al.l (j2004[ ) 
published morphological parameters, including the elliptic- 
ity, of 31 Galactic open clusters, residing basically in the 
Galactic anticentre direction. They used the 2MASS cat- 
alogue and segregated clusters from the background with 
the help of an equidensity method. We compare our results 
with the above studies in Fig. [T^l 

Our results for the two nearest clusters (e — 0.16 ± 
0.09 for the Pleiades, and e = 0.14 ± 0.09 for Praesepe) 
are in good agreement with the ellipticities derived by 
IRaboud fc Mermilliodl ^98a) and IRaboud fc Mermilliodl 
(|1998b) . This can be expected since both studies are based 
on the catalogues of the Hipparcos-Tycho family, and the 
member lists of Raboud and Mermilliod practically coin- 
cide with our membership for these clusters. Because the 
Pleiades and Praesepe are at small distances from the Sun 
and show large proper motions, one can reliably determine 
membership in these clusters up to relatively large distances 
from the cluster centres and down to stars of relatively low 
masses. Therefore, these clusters are the best candidates 
to get realistic ellipticities from observations. Indeed, the 
complete models (MC) provide e = 0.30 and 0.12 for the 
Pleiades and Praesepe, respectively, and fit the observations 
reasonably well. 
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Figure 12. Comparison of our results with data from 
the literature. The left panel compares ellipticities 
of selected clusters. The clusters from IChen et al.l 
are shown with open circles, the Pleiades 
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( Raboud fc Mermilliod 'l998a) with a red square, and 
Praesepe (jRaboud fc M ermilliod 1998b) with a blue trian- 
gle. The dotted line is the bisector. Right panel: compar- 
ison of distributions with e for 21 clusters with t > 50 
Myr within the Galactic longitude range 144°. ..250° of 
IChen et all (|2004D - filled histo gram, and our 27 clusters 
within the same longitude and age ranges (hatched his- 
togram) . The dashed (red) histogram is the distribution of 
full models (set MC) of 27 clusters. See details in the text. 



A comparison with the results of I Chen et al.l ([2003) is 
expected to be more informative since different techniques 
and different input data are used for ellipticity determina- 
tion. Unfortunately, both studies have only 13 clusters in 
common, and 12 of them do not belong to the basic sample 
since the number of their most probable members is less 
than 20 stars in our data. Nevertheless, the agreement be- 
tween the both results is reasonable: the ellipticities of 11 
of 13 clusters deviate from the bisector by less than one 
rms-error in Fig. [T2| (left panel). 

In Fig. [H (right p anel) we compare t he distributions of 
cluster ellipticities of I Chen et al.l (|2004j ) with our models 
(MC) and observations (OC). In order to avoid coordinate- 
dependent mis-interpretations, we consider only those clus- 
ters from our basic sample that are locat ed in the same are a 
of the sky as the clusters of the sample of lChen et al.l (j2004l ). 
i.e. I = 144° . . . 250°, 6 = -18° ... 20°. We also exclude all 
clusters younger than 50 Myr from th e comparison. The 
restricted samples include 21 clusters of IChen et al.l (j2004l) 
and 27 of our clusters with observed and model ellipticities. 
Whereas N-body models forecast an average ellipticity of 
e « 0.3 and a narrow spread of ellipticities in this Galactic 
direction, our observations show a much wider distribution 
enhanced at smaller ellipticities (e — 0.06. ..0.29 ). Provided 
that the peak in the distribution of IChen et al.l (|2004f ) data 
at low ellipticities (e « 0.10 — 0.15) is not random, we con- 
clude that their data suffer even more from the restricted 
area bias t han our ob servations do. The large spread in 
the data of IChen et a l. (2004) can be explained by larger 
random errors of their ellipticities. 

We briefly mention here results on th e ellipticities of 
Galac tic globular clusters. For the latter IWhite fc Shawll 
(|l987l ) determined shape parameters based on the equiden- 



sity contours method from the data of Palomar and SRC 
Sky surveys. They find that for 99 globular cluster the aver- 
age ellipticity is only 0.07±0.01, and the orientation angles 
are dist ributed randomly (bec ause the clusters are almost 
round) . IWhite fc Shawll (jl987f ) point to interstellar extinc- 
tion as a probable reason of such an unexpected result. 

5. Conclusion 

Based on a multicomponent analysis of coordinates of 
the most probable cluster members we have determined 
shape parameters of 650 Galactic open clusters. EUipticities 
and orientation angles complete the list of morphologi- 
cal and dynamical parameters (core sizes and apparent 
cluster radii, indicators of mass segregation, tidal radii) 
which were d etermined and analysed in a series of pre- 
vious papers (ISchilbach et all [20061 : iPiskunov et al.l 120071 : 
iPiskunov et al.ll2008a[) . 

We have carried out high resolution N-body simulations 
with the specially developed 0GRAPE code on the parallel 
GRAPE systems developed at ARI Heidelberg and MAO 
Kiev. The set of 27 cluster models (3 initial masses x 3 
Galactocentric distances x 3 rotation parameters) with a 
maximum number of particles of 40404 was evolved during 1 
Gyr. The cluster particles are initially distributed according 
to a Salpeter IMF in the mass range 0.08. ..8.0 tuq. For each 
particle, the decision on its cluster membership is made by 
comparing its potential and kinetic energies. 

The calculations of the apparent shape parameters of 
the modelled and observed clusters were carried out with 
the same technique. The selection of suitable models was 
based on initial cluster masses that were estimated from 
the earlier studies of our c luster sample (jPiskunov et all 
120071 : IPiskunov et al.ll2008ah . This enabled us to reahze a 
self-consistent approach for comparison of model and ob- 
served shape parameters using the data on 152 of the most 
populated Galactic clusters from our sample. 

N-body calculations show that all the models loose more 
than about 50% of their initial mass during the first Gyr 
of the evolution due to two body encounters only. A model 
cluster older than w 50 Myr keeps an oblate shape, with the 
major axis a tilted with respect to the Galactocentric ra- 
dius at an angle of «30-40 degrees. The model counterparts 
corresponding to 152 observed clusters have an ellipticity 
distribution with a peak at e ~ 0.3.. 0.4, and, on average, 
the ellipticity becomes even larger for rotating model clus- 
ters. 

However, the observed clusters show a significantly 
lower ellipticity, typically of about e = 0.2. A compari- 
son of observed cluster shapes with the detailed models 
allowed us to realize that lower ellipticities are only ap- 
parent. We explain this disagreement by a bias in observa- 
tions due to underestimated cluster sizes r2 based on the 
ASCC-2.5 data and used in the determination of cluster 
ellipticities. According to the models, the ellipticity of the 
central region of a cluster is small, it increases if outer layers 
are taken into account, and approaches asymptotically the 
largest value at 2-3 tidal radii rt . Particles outside the tidal 
radius contribute negligibly to the cluster mass, however, 
their spatial distribution is decisive for the determination 
of the cluster ellipticity. Due to a relatively low density 
of cluster members above the background in outer cluster 
regions, real clusters are usually observable only within dis- 
tances smaller than rt from the cluster centres. Though the 



apparent cluster sizes r2 based on the ASCC-2.5 are twice 
as large as previously published in the liter ature, the ratio 
r2 / rt w as found to be, on average, only 0.5 (jPiskunov et al.l 
I2008ah . 

In contrast, the deepness of the survey does not play 
such a dramatic role for the ellipticity determination. 
Although deeper surveys ofi^er a possibility to identify more 
numerous cluster members of lower masses and, conse- 
quently, to increase the random accuracy of the elliptic- 
ity determination, this alone is not sufficient to avoid the 
"restricted area bias" . 

We adopt the assumption of IWielenI (|1974 119851 ) on 
the preferential shape and orientation of cluster ellipsoids 
in the Galaxy and derive simple formulae giving a rela- 
tion between the apparent ellipticity and spatial location 
of clusters of our sample. Due to relatively high system- 
atic and random errors in observed ellipticities, however, 
we could not use the present observations for estimating 
the parameters of the reference ellipsoid. According to N- 
body calculations, the cluster ellipsoid has an axes ratio 
of a : & : c = 1.65 : 1.35 : 1, and it is tilted by an angle 
qxY ~ 30° with respect to the Galactocentric radius. 

In order to test the models with the real observations, 
the ellipticity determination for clusters must be based 
on sufficiently large areas (up to two to three tidal radii) 
around the clusters. From this point of view, equidensity 
methods are hardly useful to determine real cluster sizes, 
especially, due to a low density contrast of cluster members 
above the background at large distances from the cluster 
centres. The membership determination with kinematic and 
photometric criteria fits the problem better but requires 
surveys with higher accuracy of proper motions than it is 
in the ASCC-2.5. 
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